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some of the Physics associated with them and briefly comment on the possible applications. 
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1 Introduction 

The rigorous way to define an interacting field theory follows the Wilsonian perspective [1]. 
One starts from a conformal theory, with a given basis of operators, and assumes that this 
provides a description of the dynamics of a fixed point which approximates the behavior at 
very short distances of the field theory of interest (in the case of asymptotically free field 
theories, the conformal theory is actually free). One then adds deformations by turning on 
the couplings of relevant operators of the conformal theory (or vacuum expectation values 
for operators), which trigger a renormalization group flow to low energies (longer distances). 
The flow starts in the far ultraviolet from the original CFT (which is a fixed point of the RG 
flow) and in the deep infrared can either end into another non-trivial fixed point, or (in the 
case of confinement) into a trivial theory. The RG flow itself, and the equations governing 
it, encode all the most important dynamical properties of the field theory, and studying it 
in detail provides a very clear strategy to characterize the field theory of interest. However, 
these flows are typically difficult to study with standard field-theoretic tools, in particular 
because at some stage of the RG flow (if not all along the flow) the field theory may be 
strongly coupled. 

The AdS/CFT correspondence p] conjectures the existence of a dual description to a 
conformal field theory in terms of a string-theory background with anti-de Sitter geometry, 
and its validity is supported by countless cross-checks. It hence naturally offers a new 
opportunity to study non-trivial RG flows by making use of what today are called gauge- 
string dualities pj, thereby providing a very powerful analytical tool to deal with the study 
of RG flows in strongly coupled theories. Indeed, many examples of flows like the one 
described above have been constructed. The basic idea behind such constructions is that 
one finds solutions to the string equations of motion which have a non-trivial dependence 
on a radial direction p in the geometry, and interprets p as the renormalization scale in the 
RG flow of some dual effective theory. Provided the space asymptotes to an anti-de Sitter 
geometry at large p (in the UV), this procedure can be interpreted as providing the dual 
description of the RG flow of a fundamental field theory. 

A conceptually clean realization of these ideas is given by the flow from = 4 Super- 
Yang-Mills to a theory with minimal SUSY called A^ = 1* Yang-Mills. In the dual string 
description, the asymptotic space is AdS^ x 5*^ (with a flux for a RR five-form). This is 
deformed by the presence of non-trivial NSNS H-^ and RR F3 fields and as a result the p- 
dependence of the background is very non-trivial, and reproduces the features expected from 
the RG flow of the dual field theory (see also [5] for the dual description of the flow resulting 
from a very special deformation of A/" = 4). This statement can be tested, because the 
preservation of some amount of supersymmetry yields a degree of analytic control over the 
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coupled non-linear differential equations describing the ffow; and hence certain calculations 
can be performed on both sides of the duality, and compared. 

Another non-trivial and successful application of these ideas has been developed by Kle- 
banov and collaborators, and relies instead on the AdS^, x T^'^ geometry, where T^'^ is the 
base of the conifold. The dual conformal theory is in this case a minimally supersymmetric 
conformal field theory, described by a two-node quiver with gauge groups of equal ranks and 
with a given number of bi-fundamental matter fields [6]. A deformation is introduced in 
the string background that corresponds to an imbalance in the ranks of the gauge groups 
on the field-theory side of the correspondence. This yields an RG flow [7j that leads to 
confining field theories, in either mesonic or baryonic branches of the moduli space of the 
dual theory [HI |9], depending on the relations between the imbalance and the original rank. 
The downside of this construction, with respect to the previous case, is that the deformation 
does not switch itself off. Once a non-trivial RG flow appears, it does not reach a conformal 
point at high energies. Instead, the flow follows closely a line of fixed points described by 
quiver field theories whose rank grows without bound. Yet, by calculating correlators of some 
key-operators [10], one can see that the resulting system behaves in a way which resembles 
very closely a four- dimensional field theory. 

This particular set-up — that we will refer to as the Klebanov-Strassler (KS) system — 
gave numerous insights into the strongly coupled dynamics of four-dimensional field theories. 
These insights could find applications in Particle Physics and Cosmology, aside from giving 
non-trivial examples of calculable strongly-interacting field theories. 

The next stage in this descendent scale of 'conceptually clean' set-ups is represented by 
the so-called wrapped- brane models. Examples of these are the ones presented in [HI [12], 
and various others. The field theory here is obtained by a twisted KK compactification of a 
higher- dimensional field theory on a (curved) manifold. A weak-coupling analysis of the dual 
theory shows that it contains an infinite tower of equidistant modes (see for example |13j). 
revealing the existence of extra dimensions. Supersymmetry and R-symmetry are (partially) 
broken and the RG flow of the four- dimensional field theory never reaches a fixed point at 
high energies. These models can be thought of as dual to effective theories which require a 
UV completion — often non field-theoretic, after which a CFT is typically reached. See |14j . 
for a detailed presentation of each of the set-ups described above. A very interesting fact is 
that the deep-IR behavior of backgrounds in this class resembles very closely the backgrounds 
of models of the previous classes, while the major differences arise only at large p. Also, 
these models are comparatively simpler to construct and to study, which means that they 
are very suitable for applications. 

Another development that is relevant to this paper goes in the direction of constructing 
the duals of field theories with matter content in the fundamental representation of the gauge 
groups. Adding Nf fundamental fields in a field theory with gauge group SU{Nc) is by now a 
well-studied problem in the context of gauge/string duality. Reviews of different dynamical 
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regimes, depending on the ratio jf, their string-theoretical description and apphcations 
can be found in [15] -[17]. It is in general very difficult, within this context, to construct 
backgrounds with the nice UV-asymptotic properties of the dual to a fundamental field 
theory. 

In this paper we will construct a set of solutions that fall within the class of Klebanov- 
Strassler-like backgrounds, in the sense that the string background is very closely related 
to the KS one far in the UV and deep in the IR (i.e. the backgrounds have the same good 
features as the KS ones at the extrema of the RG flow). The novelty is given by the fact 
that the presence of Nf sources, being Nj ~ Nc, modifies the ranks of the gauge groups, and 
drastically changes the RG flow at intermediate scales. 

We will also see that such backgrounds are non-trivially related, at least at the technical 
level, to the wrapped-brane models, in the sense that they can be constructed by starting 
from one such model, modified by the addition of Nf sources, and by applying an algebraic 
procedure that amounts to a partial UV completion. Notice however that the actual back- 
grounds of interest in this paper can also be constructed directly, without ever referring to 
the wrapped systems, and that in this sense one might think of the reference to wrapped- 
brane systems as purely technical. But it is suggestive, at least from the point of view of 
effective field theories, that finding these solutions is technically much easier by starting from 
the wrapped-brane systems. 

1.1 General idea of this paper 

Using a solution generating technique presented in [18] and further developed in [T9]-[2T] 
(sometimes referred to as 'rotation') we construct a solution that generalizes that of KS [8] 
and the so-called baryonic branch of [22] by the addition of source D5/D3-branes in the 
bulk that translate into a new numerology for the field-theory quiver. The field theory, as 
we will see, makes transitions between the numerology characterizing mesonic and baryonic 
branches and shows interesting phenomena. From the view point of the string background, 
it is completely smooth and solves the equations of motion of type IIB string theory coupled 
to the action of (smeared) sources. We use techniques developed in previous works to give a 
'profile' to the sources avoiding any singular behavior. 

We analyze the field theory, whose dynamics includes two cascades: one is the familiar 
cascade of Seiberg dualities that is present in fS], while the other is associated with each 
crossing of a source D3-brane in the radial direction (corresponding to a Higgs mechanism) 
as follows from the ideas of [23]. We present various matchings between field-theory expec- 
tations and calculations with the string background. We also speculate about a possible 
interesting application to the phenomena of Tumbling in theories of Extended Technicolor. 

In practical terms, what we are doing is the following. We start from the wrapped-D5 
system, for which it is comparatively simple to find background solutions. We add fiavor 
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branes, using a profile for tlie sources sucli tliat tlie deep IR and far UV are virtually unaf- 
fected, and solve the resulting background equations. We then apply the rotation procedure, 
in such a way as to cure the ills of the wrapped system in the far UV, without affecting the 
deep IR. The result is a very novel dynamical model, which yet preserves the good features 
of older constructions. 

Because we believe the results we obtain are very interesting, in spite of the fact that at the 
technical level what we did reduces to putting together a set of well-known ideas, the work is 
presented in full detail and a fair amount of background material is also given, which makes 
the paper self-contained, if long. In a companion shorter paper [24j, we explore a particular 
realization of the ideas described above, in which the two cascades occur at different scales. 
The organization of the manuscript is the following: in Section |2l we present in detail 
the system of wrapped branes that is at the heart of our treatment. We also review the 
application of the solution generating technique of [18], ^0\. The techniques that allow to 
give a profile to our sources are described in Section O Some new solutions are displayed 
in that same section while some other new solutions to the type IIB equations with sources 
are discussed in Section HI Field-theoretical quantities are analyzed in Section |5] and some 
applications are discussed in Section |6l We close the paper with conclusions in Section [71 
Generous appendices complement the presentation. 



2 Presentation of the SUSY system 

In this section, we briefly summarize some well-established aspects of two particular four- 
dimensional supersymmetric field theories and their dual backgrounds. The two theories are 
the Klebanov-Strassler [8], [22] QFT (referred below as 'theory I') and the theory obtained 
when wrapping D5-branes on the two-cycle of the resolved conifold [12], [13] (which we 
will refer to as 'theory IF). We assume in the following that the reader is familiar with 
aspects of these two systems, and what follows is just a succinct account of their properties; 
some reviews are [H]. The theory I is a quiver with gauge group SU{n + N^) x SU(n) and 
bi-fundamental matter multiplets Ai, with i, a = 1, 2. The global symmetries arqll 

SU{2)l X SU{2)r X U{1)b x U{1)r . (2.1) 

There is also a superpotential of the form W = ^eijeaptr[AiBaAjBp\. The field theory is 
taken to be close to a strongly coupled fixed point at high energies. The dual description is 
given by the Klebanov-Strassler background [S] and its generalizations [22] • The field theory 
II has one gauge group SU{Nc) and the global symmetries are those of eq. (l2.ip — except for 
the baryonic symmetry that is not present in this system. These two theories, apparently 
so different, can be connected as discussed in [TS] and [H] via Higgsing. Indeed, giving a 

^The R-symmetry is anomalous, breaking C/(l)fl — ?• Z2Na- 
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particular (classical) baryonic VEV to the fields {Ai, Ba) and expanding around it, the field 
content and degeneracies of [13j are reproduced — see [IBj . This weakly coupled field theory 
connection has its counterpart in the type IIB solutions dual to each of the field theories. 
Indeed, it is possible to connect the dual backgrounds of field theories I and II, using U- 
duality [18]. This connection was further studied in [19], [20], [21]. In order to lay-out some 
background formalism and make manifest the connection between theories I and II at strong 
coupling, we start from the type IIB configuration describing the strong dynamics of the 
'field theory 11' (a twisted compactification of D5-branes to four dimensions). It consists of 
a metric, dilaton $ and RR-three form F^. A quite generic background of this kind can be 
compactly written using the SU{2) left-invariant one-forms 

ij I = cos ip d9 + sinip sin 9 dip , 0)2 = — sinip d9 + cosip sm9 dip , 0)3 = dip + cos 9 dip , (2.2) 
and the vielbeins 



E''^ = eidx\ EP = e^^^dp, E^ = e^^^d9 , E^ = e^-^^ sm9 dip , 

20s 9 dip 

(2.3) 



E^ = ^e^+%ui + ad9) , E^ = ^e^+%U2 - a sin 9 dip) , E^ = ^e^+^u^ + cos 9 dip) . 



In terms of these, the background and the RR three-form, in Einstein frame, read 
10 



dsl = Y.^E\ 

(2.4) 



rin2 
i=l 



where we defined 

E'^f'-i = E' A E^ A E'' A ... A E\ 



h = -2N,e~'-^^ , h = ^^-'-''{a' - 2ab + 1) , (2.5) 

/3 = N,e-'-'-\a -b), U = ^e~'-^-%' . 

The system has a radial coordinate p, on which (a, 6, (7, h., k) depend, and we have set 
a'Qs = 1. The full background is then determined by solving the equations of motion 
for the functions {a,b,^,g,h,k). A system of BPS equations is derived using this ansatz 
(see appendix of reference [25j). These non-linear and coupled first-order equations can be 
arranged in a convenient form, by rewriting the functions of the background in terms of a 
new basis of functions that decouples the equations (as explained in [26]- [2?]). We quote this 
change of basis in our Appendix \M Using these new variables, one can solve for all of them 
except for a function P(p), that is determined by a second-order equation, see [26j and our 
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Appendix |X] for details. The second-order equation reads 

P» + P'(^ + f;_|_4coth(2p,)^0, 

Q(p) = (2p coth(2p) - 1) . (2.6) 

We will refer to eg. (12. 6 p as the master equation: this is the only equation that needs solving 
in order to generate the large classes of solutions of type IIB dual to 'field theory IF in 
different vacua, deformations of the Lagrangian, VEV's, etc. 



2.1 Aspects of the SUSY solutions 

Let us start by describing some known solutions of the master equation (12. 6p . As a first 
example, the simple solution P = 2NcP gives the background of [I2], [28]. This solution and 
those with the same large-p asymptotics will not be the main focus of this paper. Other 
interesting solutions can be found semi-analytically (with asymptotic series expansions and 
numerical integration). Let us discuss the large and small-p expansions of the function P{p) 
that will be of interest in this work. For large values of the radial coordinate (describing the 
UV of the field theory II), the solution has an expansion of the form. 



e4p/3 



-i6p/3 / 18567 2781 27 ^ 2n«/3xl / 

Notice that this expansion involves two integration constants, c+ > and c_. The back- 
ground functions at large p are written for reference in Appendix El Regarding the IR 
expansion, we look for solutions with P — as p — )■ 0, in which case we find 

15 V hi ) 525 V 3/i? j'^ ' ' ^ ' 

where hi > 2Nf. is again an arbitrary constant; there is another integration constant, -P(O), 
taken to zero here, to avoid singularities. In order to match the UV expansion, it amounts to 
taking as well c_ = in (12.70 . This gives background functions that are quoted in Appendix 
\K\ Of course, there is a smooth numerical interpolation between both expansions. As one 
can imagine, there is then only one independent parameter; given a value for one of the 
constants {c+,/ii}, the requirement that the solution matches both expansions is sufficient 
to determine the value of the other. The small-p behavior of the expansion above is that 
of the exact solution P = 2NcP (but notice that in eq. (l2.8p the constant hi > 2Nc). In 
contrast, the large-p expansion P ~ e^^^^ greatly differs from the linear behavior. Let us 
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focus our attention on the dilaton field in backgrounds like those in eq. fl2.4p . For small and 
large values of the radial coordinate, we have 

^ 9/1? ^ 405/it ^ 

3Ar2e-8p/3 

1 + (1 - M + Oie'^'"'^) , (2.9) 



p4$| _ p4'I'(oo) 



^2^! p' + ■ J,. . P'QQ' = . (2.10) 



and there is a relation between $(0) and $(oo). Notice that the dilaton (that is related to 
the warp factor) asymptotes to a constant at large values of the radial coordinat^. One can 
see [19] that this type of solutions corresponds to the addition of an irrelevant operator in 
the Lagrangian of the theory of the type II field theory, needing a UV completion. 

2.1.1 A nice way of writing solutions 

We summarize here a curious way of finding an exact and analytic solution to the master 
equation ( 12. 6p . in terms of an infinite series of non-explicit integrals. The usefulness of this 
formal solution will become clear in the following sections. We rewrite the master eq. fl2.6p 
as 

sinh^(2p) / sinh2(2p)" 
We integrate this expression twice, taking already into account the asymptotics described 
above, to get 

P^-WQ^+Q r dp PQQ' -12 r dp sinh2(2p) / dp ^ = AXh^ ( dp smh\2p) 

Jo Jo J-p smh (2p) \Jq 

(2.11) 

2 4 

where we have defined for convenience the constant A = 23c+e~3, where e will be later 
identified with the deformation parameter of the deformed conifolci§. We propose a solution 
in an inverse series expansion in the dimensionless constant A: 

Pi P 9 P 

P = APi + Po + ^ + ^ + ^ + ...- (2.12) 
Plugging this expansion in (12. lip and matching the different powers of A, we get 

\ 1/3 



Pi = (^46^^" dpsinh2(2p)^ 



P0 = P_2 = ... = P-2fc = 0, (2.13) 

1 / f f p'oo' 

^-1 = -152 - P^Q^ + 2 / dpPiQQ' -A dp sinh2(2p) / dp -^^It^ 
Pi ^ Jo Jo Jp smh (2p) 



^Some readers may find this reminiscent of what happens when 'keeping the 1' in the D3-brane warp 
factor. 

■^The deformation parameter e of the conifold is such that e^^^ has dimension of length squared, the same 
as c_|- (restoring units of gsCx'). Thus A is a dimensionless parameter. 
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A recurrence relation for P-(2k+i) can be found in eq.(B.7) of [26]. This series converges 
rapidly to the numerical solutions of eg. (12. 6 p and is a very good approximation for large 
values of the radial coordinate. 

An interesting limit to take is A — ?■ oo. We will explore this in the following sections. 



2.2 Connecting to Theory I 

As explained in [19j, the solution in eqs. fl2.7p -f l2.8p . supplemented with a suitable numerical 
interpolation, describes the strong dynamics of the dual field theory II in the presence of 
a dimension-eight operator inserted in the Lagrangian (which ultimately couples the field 
theory to gravity) . This calls for a completion in the context of field theory which is achieved 
with the U-duality of |18] (a solution generating technique sometimes referred to as 'rota- 
tion'). After the U-duality described in [18] is applied, we define the new vielbein (which we 
use in the following), 

e""' = e^h~^dx\ e'' = e^^^h^dp , = e^^^h^dO , e'^ = e^'^'^'h^ sin 6 dip , 



gi = -e-+^h-^{u^ + ade), 
= ^e?+^/i^ (ws + cos 9 dip) . 



= -e 4 +^/i4 ((^2 — a sin 6* dip) 



(2.14) 



The dilaton and RR three-form are invariant under this operation. The newly generated 
metric, RR and NSNS fields are 



ds% 



We have defined 



10 

1=1 



J\2 



fli/^ L 



B2 = K 



62^ 



5 $ 



- cosa(e'^ + e^') - sma{e'^ + e^^) 
- fie''"' - /2e^'' - fsie'^' + e^^O + Me''^ + e^^S) 



-K 

-K -^dx'^ A dx^ A dx"^ A dx^ 
h 

,2$ 



h 



COS a 



cosh(2p) — a 
sinh(2p) 



sm a 



sinh(2p) 



h = 1 — K e 



2„2$ 



(2.15) 



(2.16) 
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where k is a constant that we will choose to be k = e~*'-°°'', requiring the dilaton to be 
bounded at large distance^. The rationale for this choice is to obtain a dual QFT that is 
UV-complete, see [19] for explanations. The background in eq. fl2.15p has the same form as 
the one describing the baryonic branch of the KS theory [22]. Indeed, one can check that 
the BPS equations written in [22] are exactly equivalent to our eq.f l2.6p . Another way of 
understanding the connection was explained in [20] : what relates the backgrounds in eqs. (12. 4 p 
and fl2.15p is a rescaling of the almost Kahler and complex structure forms describing the 
six- dimensional internal space. 

A closer inspection strongly suggests that the formal solution written in Section 12.1.11 is 
the solution of [22]. Let us study the limit A — oo of the special solution presented in Section 
12.1.11 In that case the solution reduces to 

P ~ APi = (4e^y" sinh2(2p)^ ' = A 2-^ei (sinh(4p) - 4p)i/^ (2.17) 

For the different background functions the limit results in 

2k pi 1 



2p 4<E.-4$o = 2sinh^(2p) 

sinh(2p)' \^P^P[ ^ ^ ^ 



/ 2sinh^(2p) , V^" sinh(2p) (2PiP{P^i + PfPi, - P-Q^) , 

The choice of k = e~**^°°-* implies that = 2e''A^e~^*°, and the logic for this choice is 
to get a warp factor that vanishes at large p. Indeed, plugging this value in eq. fl2.18p . and 
using the explicit expression of Pi given in eq. fl2.17p . we see that the two first terms in the 
expansion for h cancel out and 

y P^^/3Ff ^ ' \' ^ ' ^ ' 

This in turn implies that we are switching off a dimension eight operator in the dual QFT, 
obtaining a proper UV-complete field theory. A detailed discussion of this choice, showing 
that the decoupled operator has dimension eight, is in [19]. Regarding /i2 we explicitly have, 

= ^^N^ r ,,(W-4x)(2xcoth(2^-l) _ ^^^^^ 



p8/3- 



sinh^(2x)(sinh(4a;) - 4a;)2/3 



*This is the restriction on the type of solutions considered that rules out the exact solution P = 2NcP- 



10 



which is the Klebanov-Strassler [8] warp factor h2 = h^s- Finally, rescaling the Minkowski 
coordinates — )■ x^A"^/^ the metric is independent of the parameter A, and using that 

Pi_ 8sinh^(2p) 

Pi 3(sinh(4p) - 4p) ' ^ ■ ^ 

the internal space metric is the deformed conifold. This is precisely the Klebanov-Strassler 
solution [8j that is obtained as the limit A — )■ oo of the formal solution in Section 12.1.11 (the 
latter is the full baryonic branch solution [22j). 

2.3 Adding sources 

We consider now the addition of sources in these backgrounds and the effect of the rotation on 
them. The logic followed is the one described in [25], [16]. We add a bunch of Nf D5-branes 
in backgrounds of the form given in eg. (12.41) . These D5-branes are added as explained in [25] . 
They are in principle localized in the compact space described by the angles [9,ip,9,ip]. One 
can run the machinery described in [20] and [29], to write a BPS system of nonlinear coupled 
partial differential equations (with derivatives on the radial coordinate and the angles where 
the sources are localized). The difficulty of the problem prompted the authors of [30], [25] . 
to smear the source^. In this smeared, SU{2)i x SU{2)r invariant situation one can run a 
very similar program to the one described above: generalize the change of variables discussed 
in Appendix |A] and write a new master equation for the function P{p), 



0, 

Q(p) = ^ — ^(2pcoth(2p) - 1) . (2.22) 



2Nr - Nf 



See the paper 126] for technical details. In this case, solutions are known either semi- 
analytically (as series expansions plus a smooth numerical interpolation) and also gener- 
alizing the formalism of Section 12.1.11 A classification of the solutions according to the 
asymptotic behavior of -P(p), together with expansions of all other background functions 
in the presence of sources was given in [26]. We will not quote the results here, but refer 
the reader to [26] and [Sn] for details. A feature of all these solutions is the presence of 
singularities for small p. In other words, one finds a divergent Ricci scalar as i? ~ A 
physical reason of the existence of the singularity (namely a very large density of 'crossing' 
source-branes at the point p = 0) was given in [16], among other papers. More in detail, 
the solutions in the presence of the D5 sources are formally like those in eq. fl2.4p . with the 



""Equivalently, either consider the Fourier zero mode of the solution in the compact internal space, or 
consider a situation where the global symmetries of the field theory without fiavors/sources — eq. p.ip — 
are respected after the addition of sources. See [16] for discussions on different aspects of the smearing. 
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difference that the functions fi, f2, fa, fi change to reflect the presence of sources {dF^ is 
nonzero). Indeed, we now have, 

/, = -2iV, e-^-^^ , f, = ^e-'-'\a' - 2a6 + 1 - ^) , 

h = N,e-'-'-^{a -b), U = ^e-'-'-^ b' . (2.23) 

The other functions in the background also reflect the presence of the Nf sources. After the 
rotation procedure is apphed we will have a background that formally is the one in eq. fl2.15p . 
but with the functions $, h, g, k, a, b that now have contributions coming from the sources. 
There are two interesting points to discuss. First, as explained in detail in ^20j, the presence 
of the NSNS B2 field will induce D3-brane charge on the D5 sources after the rotation. The 
second point is that due to the presence of D5 sources the function P{p) gets a contribution 
at large values of the radial coordinate, see [26] and compare with eq. fl2.7p . 

P(p) = e^^/^ (^c+ + ^e-^p/^ + 0(e-«^/3) j _ ^2.24) 

Carefully following the details — see [20] — one finds that this new term impacts the large- 
radius asymptotics of the warp factor for the rotated background as 

h = Nfe-^P/' + e-8p/3 ^ 0(e-8p/3) (2.25) 

2c-|- 

Then, asymptotically, the cascading behavior of the KS system, represented above by the 
g-8p/3 ^Qj-m present also in h2 of eg. (12. 20 p . is overcome by the presence of the D3 and D5 
sources. This deviates the system from the 'near AdS^ UV. 



2.4 About the dual field theory 

We comment briefly on the field theory dual to the (singular) geometry after the rotation 
with sources. This was proposed in [20] to be the quiver 

SU{n + N^ + Uf) X SU{n + nj + ^) , (2.26) 

{uf is the number of D3-branes induced by the sources, more precisions to be given in the 
following sections) with superpotential W ~ ABAB for the fields Ai,Ba transforming as 
bi-fundamentals of the gauge groups, in the same fashion as indicated around eq. fl2.ip . It 
was also proposed that the field theory was on a mesonic branch and that as discussed in [9] 
the baryonic symmetry U{1)b was unbroken. Computations of the c-function, beta functions 
and anomalies supported this proposal. The presence of the small-p singular behavior of the 
background (inherited from the singularity in the unrotated background), implies that we 
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should not trust the IR dynamics as read from the geometry. At best, the solution should be 
thought of as a generalization of the Klebanov-Tseytlin geometry [7] in the presence of extra 
matter/sources. This point was made concrete by calculating Wilson loops that displayed a 
non-physical behavior for large separations [3T| . 

It was clearly stated in [20] that a resolution of the singularity was very desirable as 
it would allow to trust the low energy dynamics of the field theory read from the type IIB 
background with sources. Here, we also point to the need of having a better understanding of 
the high-energy (and strongly coupled) dynamics. As pointed out above, this dynamics seems 
to be driven by higher-dimensional operators reflected in the term Nf e'^^l"^ in eg. ( 12.25^ . We 
will comment and resolve both issues in the following sections. 



3 Solving the problems: sources with a profile 



Let us briefly remind the reader about the recent results of [22] • The paper is written in the 
context of the backgrounds of eg. (12.41) for (iVc, iVj) D5 color and flavor branes. The authors 
found that a way of getting rid of the small-radius singularity was to avoid the 'very high 
density' superposition of sources (as described above and also pointed out in [16], [33]). For 
this, the authors of j32] developed a formal way of distributing the sources such that not all 
the N f branes reach the point p = 0. They constructed a distribution of sources (flavors) 
with a particular profile that vanishes near the origin and stabilizes to a constant at large 
radius as indicated in the Figure [TJ 



D5-brane 




Figure 1: We plot an example of a possible D5-brane configuration: a total of N ^ D5s extend 
in the radial coordinate from some minimum p, that can lie between p^, and p**, to infinity. 
We chose a distribution of minimal radial distances of the D5s that is homogeneous in p, 
as indicated by the dotted obligue line. This configuration of D5s gives rise to the plotted 
^(p), which is asymptotically constant {S{p) — ?■ 1). 



The distribution S'(p) is obtained rigorously solving the kappa symmetry eguations for the 
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branes and distributing them following a well-defined procedure [32], [35 • This distribution 
can be interpreted in two ways: the first, as the addition of radius-energy-dependent masses 
for the flavors in the dual field theory — the profile of masses is somewhat related to Nf S{p). 
The second way of interpreting the distribution of sources, uses that the field theory II couples 
to flavors according to the superpotential (see [25], [27] ) 

W^Q^^kQ, (3.1) 

where [Q) Q are the (anti) quark superfields introduced by the sources and is a generic 
massive chiral multiplet in the field theory II. The profile is understood as an energy- 
dependent VEV for the field that effectively generates an energy-dependent mass term 
for the flavor multiplets. While the outcome in both cases is an energy- dependent 'distri- 
bution of mass', the second interpretation implements it without breaking the R-symmetry, 
hence making it our choice. Let us describe how this distribution of sources affects the string 
background. 

Indeed, one writes a background like that in eg. (12.41) . where the functions /i, . . . , /4 now 
reflect the presence of smeared sources with a profile S{p) B 



A = -2iV,e-'=-^^ , h = ^e-'-"^ (a^ - 2a6 + 1 - ^5(p) 

^' = ^^' ^' = Y' l,^+2iV.cosh(2p)J 



(3.2) 



Following [32] , one can write the BPS equations for the system with the sources added with 
a profile. Here again, the formalism of [26] can be generalized, in this case for a distribution 
Nf{p) = NfS{p). Indeed, all functions can be immediately obtained once we know the 
functions P, Q. Following the treatment in [32] we learn that 



wo N i" , 2Nc- Nf Six) 



where an integration constant has been set to zero to avoid IR singularities. Moreover, as 
shown in [32], one can find a master equation, 

which reduces to eq. fl2.22p when S{p) = 1 and to eq. fl2.6p for Nf = S{p) = . 

Once the function S{p) is known (from a microscopic description of the smearing of 
sources), one can get Q{p) from (13. 3p and solve the second-order master equation (13. 4p 



^Compare this to eq. (|2.23p to notice that the effect of the profile for the sources is not just the change 
Nf NfS{p). See ^ for details. 
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for P(p). The rest of the functions have expressions that generahze those in our Appendix 
Rl as we read from 



-,2h 



-.2k 



4 P coth(2p) - Q ' 
P' + Nj S{p) 



=29 



P coth(2p) - Q , 
P 



(3.5) 



The dilaton is, 



^4$-4$, 







Pcosh(2p) -gsinh(2p) 
2sinh(2p)2 



and the function b{p) has the form 

2p Nf 



b{p) 



{P^ -Q^){P' + NfS) 
S{p) , 2 



sinh(2p) 2N^ Lcosh(2p) sinh(2p) 



dx tanh(2x) S{x) 



(3.6) 



(3.7) 



The problem just boils down to finding P{p), given S{p). We follow the treatment in 
[32] and |3l] — where quite generic distributions S{p) were constructed — to get particular 
expressions for the function S{p) that make the expression of Q{p) analytic and the resolution 
of the master equation simple. Performing the rotation is immediate, and we will quote 
below the expression for the full rotated background. As we discussed, the rotation induced 
D3-brane charges on the D5 sources. This is represented in Figure El 

Nf S{p) 



Nf 









- D5-brane 






X 


D3-brane charge 


;a 


_x 
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Figure 2: After the rotation, D3-brane charge is induced on the D5s. As argued in Section 
13.11 this charge grows exponentially with the radial distance. An asymptotically constant 
profile S{p) will generate an infinite amount D3-brane charge. 



3.1 The two kinds of solutions 

Given a distribution function for the D5 sources S{p), we can actually construct two kinds 
of solutions: those that are exact and analytic and those that are semi-analytic (in the sense 
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of knowing either their asymptotic power series expansions or an exact expression in terms 
of an infinite number of integrals that cannot be evaluated exactly). 



• Exact solutions: they are obtained in the double-limit A — ?■ oo and Nf — )■ 0, keeping 
XNf fixed. We then define z/ = 2-'^/^e^/^XNf = c+Nj, where c+ was the dimensionful 
constant in equation (12. 7p . The solutions are very similar to the KS background in 
the sense that functions are written up to an integral, the dilaton is constant and 
F3 = *6-f^3- In these solutions there are D3 sources, but no D5 sources. 

• Semi- analytic solutions: the parameters A, Nf in the expansion are both fixed and 
nonzero. These solutions are written in terms of an infinite sum of integrals. The 
treatment follows what is written in Appendix [Bl For all practical purposes, this is 
equivalent to writing asymptotic series for large and small values of p and a smooth 
numerical interpolation. In these solutions there are D3 and D5 sources. 

We will give details on the semi-analytic solutions in Section HI In the case of exact solutions, 
the configuration in string frame reads (details of the derivation are given in Appendix [B]) 

ds^ = h-'/^dxl, + h'/^dsl, 

dsl = Pi (^{dp' + ]{CJs + cos ed^f) + ^^^y p) ^ ^.^2 ^^^2) ^ 

V 2x1 4 4 

+ S^^MM ((^^ + adey + {C02 -asinedip)') 

B2 = _:^ ^/^^"^^^(^/^) — 1 (cosh(2p) (sin 6 d6 A d(p - ui A CJ2) - sin 9 d(x> A Coi - dO A UJ2) 
4 smh(2p) 

Hs = dB2, F3 = HH3 , 

F5 = -(1 + *)dph~^dx° A dx^ A dx^ A dx^ A dp , 
It ^ r, S{x) 

h = h2 + I dx "T^r , 

e^/^ Jp (-4x + sinh(4x))^/^ 



cosh(2p) 

(3.8) 

where Pi is written in eq. (l2.13p . and /i2 = hxs is the Klebanov-Strassler [8] warp factor, as 
written in eq. (l2.20p . The dilaton is a free constant and we take e* = 6*'-°°-' = 1. Let us 
stress that when one takes S* = in eq. (l3.8p . one obtains precisely the Klebanov-Strassler 
background. 

Serendipitously, we write the warp factor as 

^ ^ 4 g(x)(-4a; + sinh(43;))^/^ 

e ' J p [—Ax + smh(4x)) ' 
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The function 

f oo 2 



GiP)= I (3.10) 
(—4a; + sinn(4x)) ' 



p 



indicates that the solution above is just the solution to the Laplace equation (for the function 
h) in the presence of fluxes F^.H'i.F^ and D3 sources (see eq.(lOl) in [8]). So, we have the 
KS geometry being deformed by a distribution of source D3-branes 

~ z/ S{p) (sinh(4p) - 4p)^/^ , (3.11) 

that are supersymmetric when placed on the deformed conifold. Hence, this exact and 
analytical solution could have been written without going over all this effort, just assuming 
the strange distribution of D3 sources in eq. fl3.11l) . On the other hand, the solutions with 
A, Nf finite -presented in the next section — are new. Since the asymptotics of both kinds 
of solutions are quite similar, we can draw generic lessons just by studying the exact analytic 
solutions in eq. (13.81) . 



3.1.1 Generic features of the solutions 

Let us analyze a couple of interesting points. According to what we wrote in eq. (l3.1ip . after 
the rotation, the induced source D3-branes are distributed in such a way that they pile up 
exponentially towards large values of the radial coordinate as n/ ~ e^''/^. This is precisely 
what produces the solutions' departure from the four dimensional behavior of the cascade. 
This large pile-up of D3-branes dominates the UV dynamics and is equivalent to the insertion 
of a dimension-six operator into the Lagrangian, as we can see by expanding the warp factor 
h in eq. (l3.8p giving as a result the one in eq. (l2.25p . We will analyze later how to get rid of 
this undesirable behavior. 

The second interesting point (as anticipated above) is that for a suitably chosen distri- 
bution S'(p), we can make the number of sources decrease towards p = 0. In this way, one 
avoids the presence of the small-p singularity. We can follow the papers [32]- [3l] to choose 
an adequate S{p). The backgrounds constructed present no pathology at p ~ 0. Hence the 
low-energy strong dynamics of the field theory I can be calculated using the backgrounds 
above and is fully trustable. We still need to do something about the UV behavior, that 
takes us away from the cascade behavior. Let us make some comments on the field theory 
dual to our backgrounds that will illuminate the way to resolving our problem in the large 
radius region. 



3.2 Field theory comments 

In [9] it is explained that one may think the moduli space for the quiver 

SU{N^+p) X SU{p) , (3.12) 
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as {p — INc) D3-branes free to move on the deformed conifold with ((/ + l)Nc)D5 together 
with INc anti-D5-branes forming a bound state at the threshold. The ways of distributing 
the {p — INc) D3-branes give rise to a symmetric product of deformed conifolds as moduh 
space. An important point made by |9] is that if we start from a quiver SU{{k + l)Nc + p) x 
SU {kNc +p) — in the mesonic branch of the KS field theory — and Higgs down the group 
repeatedly, we will arrive at the quiver SU{{k + l)iVc) x SU{kNc) x ?7(1)^. This quiver has 
the numerology to fit the baryonic branch of the KS field theory. The difference is that the 
U{X)baryonic IS gauged as it mixes with the [/(l)^. 

The authors of [9] study carefully the quantum dynamics of the quiver in the various 
ranges of p, taking into consideration the presence of the tree-level and the quantum-induced 
superpotential. They show that the deformation of the conifold is always proportional to 
the scale at which the theory with larger gauge group goes strongly coupled. Let us now 
apply this information to our case with sources. 

Let us consider for the following subsections the case in which the function S{p) vanishes 
for p < and stabilizes to S{p) — )■ 1 for large values of the the radial coordinate. Also, 
let us focus our attention on backgrounds like the one in eq. fl3.8p . though the lessons will be 
valid also for the semi-analytic solutions that we present in Section HI 



3.2.1 The theory at high energies 

We study the UV of the field-theory dual to the background in eq.f l3.8p . We interpret this 
solution combining the results of the papers [9] and [20]. We propose that the dual quiver is 
of the form 

SU{n + nf + Nc)yiSU{n + nf), n = kNc. (3.13) 

The field theory is in the mesonic branch for p > p^,, where the function S{p) ~ 1. There are 
two competing processes, as explained in [20], inspired in early ideas presented in [23]. The 
usual cascade, represented by the hxs in the warp factor of eq.f l3.9p and a Higgsing process 
represented by the term proportional to u in eq. fl3.9p . In the new radial coordinate r ~ e^^^'^, 

hlp^oc ^ , (3.14) 

we see clearly the superposition of the cascade and the Higgsing. The presence of an expo- 
nentially increasing number of source D3s effectively behaves as the insertion of an irrelevant 
operator of dimension six, deforming the UV dynamics out of the near-conformal KS dy- 
namics. This indicates the need for a UV completion. 
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3.2.2 The theory at low energies 



Flowing towards the IR, and close to p ~ the Higgsing and cascading lower down the 
ranks of the groups, so eventually we will reach the numerology of the baryonic branch 

SU{nf + {k + 1)N^) X SU{nf + kN^) 

SUijif + {k + l)N^ - 1) X SUinf + kN^ - 1) x f/(l) ^ (3.15) 
^ ... SU{{k + l)Nc) X SU{kNc) X f/(l)"/ . 

The main difference with respect to the usual baryonic branch as discussed in detail in [9] is 
the fact that the baryon symmetry is in this case gauged. This happens because the baryonic 
symmetry mixes with a diagonal combination of the [/(l)'s that appear when Higgsing. This 
symmetry being gauged impacts the low-energy phenomenology. For example, we can check 
in the case of the background in eq.f l3.8p . that the massless excitation described in [35] is not 
a solution to the equations of motion. This implies that the theory is on a mesonic branch. 

For the rest, with a suitable choice of S{p) vanishing conveniently fast at p = 0, the 
phenomenology of the low-energy field theory is very similar to that of the Klebanov-Strassler 
theory, as both backgrounds are quite similar. The computation of various IR quantities will 
give qualitatively similar results to those computed with the backgrounds of [8] and [22j. 
There will be numerical differences depending on the profile S{p). We will provide examples 
of this in the following sections. 



3.2.3 How to improve the UV behavior: a phenomenological approach 

We would like to stop the growing number of D3 sources responsible for the pile-up of 
Higgsing processes and ultimately for the deviation from the nearly conformal behavior 
as illustrated in the discussion around eq. fl3.14p . A way to do so would be to choose a 
profile for the sources that somewhat recovers the cascade behavior. In the case of the 
background in eq. fl3.8p this is an easy task because, as we are distributing the D3 sources 
on the deformed conifold, the solution preserves supersymmetry for any S{p). This suggests 
that we should choose — at least at a phenomenological level — a distribution function 
behaving as S{p) ~ e"^'^/^ for large values of p. A profile that does the job of fixing all the 
IR and UV problems is inspired in the work of [M] and [32], 

^(p) = tanh(2p) (3.16) 

In Figure [3] we have plotted the profile function (I3.16p . By using this distribution functioiill, 
one can compute the new warp factor using eq. (13. 9 p . We will then have a smooth IR geometry 
and a cascade behavior in the far UV. 

^Notice that one is taking p, = in the example profile p. 161) . It would be trivial to introduce the scale 
p* by writing something like S = Q{p — p*) tanh(2p — 2p*)'*e^'*(''^''*)/'^. 
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In principle, one could of course try different profiles that make the distribution of D3- 
brane sources either constant or vanishing for large values of p. However, if we want this 
distribution of sources to have a positive mass density everywhere (and vanishing for p — )■ cxd), 
it seems that the only possibility is to have profiles decreasing exactly as e"^''/^. See Appendix 
Ofor details. 

S(p) 



0.2 



0.1 




Figure 3: Plot of the profile function f l3.16p . The function S{p) starts growing for small p, 
reaches a maximum, and then decreases exponentially for large values of the holographic 
coordinate. 



Another way to motivate why the profile S{p) should decrease towards the UVu is to 
consider the Klebanov-Strassler [S] geometry with N^. D5-branes and add to it IN^. extra D3- 
branes. The field theory is then in the 1-th mesonic branch [9] of the theory I with groups 
SU{{k + l)Nc) X SU{{k + / + l)Nc). The relation between the deformation parameters e/ and 
eo of the conifold is given in P, pS] , 

27ri 2tvINc 

ei = ele^ =ele^ . (3.17) 

So, we see that 

iV^Moge; -/A^.logeo. (3.18) 

This implies (for similar deformations e/ ~ eo), that IN^ ~ N^. That is, the number of source 
D3-branes that we called nf should be of the same order as the number of flux D3-branes 
that we called n in our quiver of eq. fl3.13p . Obviously, the situations we discussed above, 
where the source function S{p) — )■ 1 at large values of the radial coordinate, generate as 
explained in eq.f l3.1ip an exponentially large number of D3 sources Uf. Without a damping 
factor like the one proposed in eq.f l3.16p or similar there is some feature of the field theory 
that is not captured by the solution. So, we will choose functions describing the profile of D5 

^This argument was suggested by Anatoly Dymarsky. 
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sources (this will directly affect also the profile of D3 sources) of the form given in eq. fl3.16p . 
These profiles decrease fast enough in the asymptotic regions, so that the IR singularity 
generated by a large density of sources and the UV appearance of an irrelevant operator are 
both avoided. 

3.2.4 Further comments on profiles and Goldstones 

To avoid confusion with the different uses and meanings of the function S{p), some words 
are in order. The function S{p) originated in the supergravity background dual to the 'field 
theory IF with D5-brane sources |32]. Its meaning there was that of a profile for the D5- 
brane charge present in the background. Coming back to Figure [H we are placing several 
D5-branes on the geometry, each of them reaching a minimal radial distance. These minimal 
distances are distributed on a shell. An electrostatic analogue of this situation is given by the 
electric field created by a hollow cavity with a thick charged shell. In this case S{p) would 
count the effective radial charge. Inside the cavity S{p) = as there is no electric field. As 
we cross the shell, S{p) increases, and stabilizes to the total charge away from the shell. The 
electric field outside the shell, will be that of a point charge. Had the shell null thickness, 
the electric field would display a jump. This is what the S{p) in Figure [1] represents. Away 
from the shell, S{p oo) ~ 1 in coincidence with the S{p) generated by 'massless' sources. 

This interpretation of S{p) as a profile for the D5-brane charge is however lost after the 
rotation. We gave an interpretation of S{p) for the exact solution of eq. fl3.8p . As explained 
after eq.f l3.10p . in this new solution S{p) accounts for how the D3-branes are distributed 
in the geometry. So, for this exact-analytic solution, S{p) represents a distribution of D3- 
branes. More precisely, the distribution is given by S{p) (sinh(4p) — 4:p)^^^. 

A curious fact is that the 'physically sensible' S{p) for the unrotated background (which 
asymptotes as ~ 1), is pathological for the rotated configuration as the number of D3- 
branes grows exponentially with p, inducing an irrelevant deformation. Conversely, a physi- 
cally sensible S{p) in the rotated background (vanishing at infinity as e~'^^^P) does not have 
— at present — a microscopic interpretation in terms of the usual flavor D5-branes in the 
unrotated solution. A clarification of this point remains for future work. 

About Goldstone bosons: Let us briefly discuss the situation with Goldstone bosons. De- 
tailed calculations needed to sort out the issue in full are left for future worlj^. 

First, let us state the question. Consider a background like the one in eq. (13.80 . We can 
choose any S{p) with a similar behavior to that in eq. (l3.16p . without affecting the argument. 
In practical terms, we are making a transition from a far UV with similar dynamics to KS 
and the baryonic branch, to an intermediate region, with the numerology of the mesonic 
branch, to a deep IR which has the gauge-group ranks of the baryonic branch of [8]. If 
one then calculates the solution for the massless mode in [35], the proposed ansatz of |35j 

^Many thanks to Anatoly Dymarsky for crucial comments reflected in these following paragraphs. 
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does not verify the linearized equations of motion. Indeed, despite the fact that F3 and B2 
in eq.f l3.8p are the same as in the KS background, the five- form flux F5 is modified by the 
presence of sources (as seen from the fact that h 7^ hxs)- Then, in particular, the equation 
of motion for would be different from the standard KS one, including a term involving 
dF^ which is non-zero because of the sources. 

In line with this, as we emphasized above, transitions between mesonic and baryonic 
branches should produce gauged baryonic symmetries p], and as such the Goldstone of the 
baryonic f/(l) should not be part of the physical spectrum (it should be Higgsed away into 
the resulting massive gauge boson). 

Nevertheless, one may speculate that a modification of the ansatz in |35] may reveal the 
existence of a massless mode. If we were to find a Goldstone boson on the gravity side, how 
could we understand it in terms of the field-theory dual? In the combination of Higgsing and 
cascading, there is a complicated pattern of global and gauge symmetries, which are non- 
trivially mixed with one another, and which undergo partial symmetry breaking. Vacuum- 
alignment arguments suggest that, if possible (i.e. in the absence of non-trivial obstructions 
within the Lie groups and cosets describing the symmetry of the system), the vacuum of the 
theory will align itself with the directions in the internal space of the global symmetry. In 
doing so, the vacuum will prefer to spontaneously break part of the symmetry which is not 
gauged. As we wrote, a detailed calculation would be needed in order to establish whether 
or not there is a normalizable massless mode in the spectrum, and we leave this for a future 
study. 

We move now to the study of the semi-analytic solutions described in the beginning of 
Section 13.11 

4 Details on the new solutions 

In this section we give details for the semi-analytic solutions. We will concentrate on the 
configuration after the rotation. We consider the case in which the parameters A^/,A are 
free. While the solutions could be written in terms of an infinite series of integrals (that 
cannot be performed exactly), we find it more useful to give series expansions close to 
p = and p — )■ cxD and a numerical interpolation between both regimes. After solving the 
master equation (13. 4p . we will use the expressions in eqs. (l3.5p - (l3.7p to quote the asymptotic 
expansions for the functions appearing in the background. As anticipated in eq. (l3.16p we 
will consider (this is just a possible example) the function S{p) = tanh''(2p)e~^''/'^. 

We should clearly state that we are not obtaining the profile above in a rigorous fashion. 
Since the distribution of D3-branes on the deformed conifold can be chosen (they are SUSY), 
the solutions with A — ?■ oo,A^j — t- of eq. (l3.8p are correct. The ones we present in this 
section, lack a rigorous proof of existence (we are choosing a profile of D5-branes that we did 
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not prove to be derived from a kappa-symmetric embedding). The healthy behavior of the 
backgrounds we found (together with the positivity of Too for the sources) suggests that the 
brane embedding can be rigorously derived. Let us start by studying the large-p expansion 
of the solutions to the master equation ( 13. 4p given the source term in eq.( 13.16p . 



4.1 Large-/? expansions 

First of all we define 

poo 

S^= / rfp tanh'(2p)5(p) , (4.1) 
Jo 

so that we can write for the asymptotic behavior of the function Q{p): 

Q = N,{2p - 1) - NfS^ + \Nje-^f'^ + (4iV,p - 2iV^^oo)e-^'' + 0{e-^''p''') . (4.2) 
The function P{p) is given at large p by 



Pip) = (c, + ^e-"/^ (4p= - 4p(l + ^SJ + H + 5£i|L + 2^ 

+ ^ {2N,N,p^ - 2pHcl - N,N, + iVjSoo)) + ©(e-'^'V"; 

(4.3) 

Using eqs. fl3.5p - fl3.7l) . we get for the background functions 

h =3^^ (iV2(8p - 1) + 2c+Nf - 4N,NfS^) - -Le-^P[2N,{p - 1) - NjSj) + 



(l6iVf p2 - 16iV,(iV, + NfS^) p + 



,2H _^ 4P/3 , iV.(2p - 1) - iV^goo 

4 + 4 + 16c+ 

+ 13Ar2 + Qc+Nf + SA^eA^/^oo + ^N]Sl}j + 0(6"^"/=^) , 



=29 ^, AT (On - ^\ - AT p-ip/'i 

2fc p-4p/3 

^ =ye^''/3 - ^(4Ar,p - SAT, - 2NjS^f+ 

^ ^-8,/3 (44 + NfS^){^p - 3) - Ar,Ar^(8p^ + 2p - 3) ^ 

36c+ 

c+ 
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2c^ 



/ Qp-8p/3 

g4* ^g4$(oo) _ ^^2(g^ _ ^) _^ 2c+iV; - AN^NfS^] + 

+ ^e-^^(2iV,(p - 1) - NfS^)) + 0(e-i^''/=') , (4.4) 

where e^*(°°) = Notice that one effect of the sources — with this given profile S{p) 

— is to add to the D3-brane charge in the warp factor. The sources do not interfere with 
the cascade, that is now driven by the bulk D3-branes only, associated with the number n 
in our quiver. Let us discus the IR expansion next. 

4.2 Small p expansions 

With the same S{p) as in (13.161) . we can see that close to p = 0, the effect of the sources 
vanish as 

64 640 

S{p) = 16p^ - ^p^ - + O(p^) . (4.5) 

It is clear that the deep IR asymptotics will suffer little modifications from that of the case 
with Nf = S{p) = 0. In order to keep track of the IR modifications, we just need to look at 
the terms with a factor Nf and compare the expansions with those in eg. ( 1A.6P and eqs.l lA.Tl) . 
The functions Q{p),P{p) read 

Pip) = hp + 4^Y^p' + Y^(3/^f - ^hlN'^ - 32iV,^ - 450hlNf)p' + 0{p') , 



while the functions in the background have expressions given b; 

„2$(0)-2*(oo) o,/o/\r2 
h=l- - ^ ' g2#(0)-2#(oo) 2 

2^2 9hl 

_ 2v|^2*{o)-2*(oo) f^^N^ ^ 224^ + 135^ ) p' + 0{p') 
135 y ii^ hi hi 



2 45 V h 
16 

+ (267/i^ + ISUhjN^ - rOOhiN^ - 328N^ - 15hl{77N, + 45Nf)) p^ + C(p^) 
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Notice that for the function a{p) above, the first correction of the sources appears at order 
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15 V ' ^ hi 



+ — - — — — + ""^'^ 



f:^hi, {hi- AND p^ ^ {Qhf - ShjN^ - - 270hlNf) p^ 
4 8 lO/ii 315/1? ^ 



4<I> 



4^(0) f, 64iV,V^ 16(-120iVf /^^ + 992iV^^ + A05hlNf)p' \ 

V 9/1? 405/if ^ V ' 

l+(-2 + ^] + 2 (75/^; - 232/^?iV, + 160h,N^^ + 64iV3) ^ 

\ 3/ii y 45/if 



& = 1 - ^P^ + (| - ^)P^ + 0{p') , (4.7) 

where 6^*'-°^ = ^w^. Using these expressions, one can find a smooth numerical interpolation 
between the two asymptotics and compute observables with them. In Figure we plot the 
functions P{p) and Q{p) for S{p) = tanh''(2p)e~'^''/^. Using P{p),Q{p) and S{p) all other 
functions in the background can be obtained. Studying these functions, we do not observe 
any pathology or non-physical behavior. Hence, we consider the choice of S{p) in this section 
to be good. 

5 Further analysis of field-theory quantities 

In this section, we continue the study of the background found in Section HI We compute 
some quantities in the proposed dual field theory. The purpose of these calculations is to 
investigate if the profile chosen for S{p) produces any undesirable feature that might suggest 
non-physicality of our arbitrary choice. We start with the central charge in the dual QFT. 

5.1 Central Charge 

We calculate the central charge c, following the prescription developed in [37], [38]. Start 
with a metric (in string frame) of the form 

ds^^ = a dx\ ^ + a 13 dp^ + Qij dy'^dy^ , (5.1) 

and following [38] we define Vint as the volume of the internal space, and H as 

H = e-'''Viy. (5.2) 

We have that the central charge (for c/ = 3) is 
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Figure 4: For the case Nc = Nj = 2 and /ii = 5, we show the functions defining the 
semi- analytic rotated background. The choice of S{p) is the one specified in the text (eq. 

(M)). 

In particular, for a background like in eq.f l2.15p . we have 



a 



^2g2<I>+2?i+2g+4fc (•54^ 



8(9pl0g[V^e2*+2'^+29+fc])3 ' 

We consider the semi-analytic backgrounds of Section|l]and plot the central charge of eq. fl5.4p 
in Figure |5l We observe a monotonic behavior. This quantity does not show any sign of 
unphysical behavior for the chosen profilJ"]. 

For backgrounds like those in eq. (12. 4p — before the rotation — we should take k = and 
/z. = 1 in the formulas above. 



^^One could have arbitrarily chosen a different profile, for example S ^ tanh''(2p)e ''"^ . In this case 
we would have observed a non-monotonic c-function. 
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Figure 5: Plots of the function c(p) for Nc = Nf = 2 and hi = 5. We considered the 
backgrounds of Section HI 



5.2 Quantities defining the gauge groups 



In specifying our quiver, for example around eg. (13. 131) . we have introduced various quantities: 
n, that is the number of bulk (flux) D3-branes; Uf, that is the number of source D3-branes; 
Nc and Nj, the numbers of bulk and source D5-branes. Here, we give a geometric definition 
for these quantities, forgetting about numerical factors. 



n + nf= I F^= B2AF3 
'25 Jt,s 



n 



Nf=0 



(5.5) 



Nf=0 



where S5 is the five-dimensional manifold of coordinates {6, (p, 6, (p, ip). Using the expressions 
in eq. (l2.15p . we can calculate F^\^ to be 



cosh(2p) — a 
sinh(2p) 



-2Nre 



2h 



Nc 



e^^a^-2ab + l] 



N 



f . 



AN, 



sinh(2p) 



e'^ia - b) 



2 ' '2 
sin 9 sin 9d9 A dip A d9 A d(p A dtp , 



5 - 



(5.6) 
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that translates into 



2$ 



cosh(2p) — a 
sinh(2p) 



+ 



2h, 



n 



sinh(2p) 

{n + nf)\Nf=o, nj 



2iV,e2^ - ^e^9^a^ _ 2ah + 1) + ^e^^S 



a-h) 



(5.7) 



{n + Uf) - (n + nf)\Nf=o ■ 



While it is good to have the exact expressions above, for many purposes the large and small 
p expansions are illuminating. For low energies in the dual QFT, we have 



n (X NIp^ + 



( ^hl 



Uf (X Nfhi 



which makes it clear that deep in the IR the Higgssing (represented by n/) has stopped. 
The geometry is very similar to that of [22] and so is the dual QFT. To describe the regime 
of high energies in the field theory, we change to a radial coordinate r = e^''^^ as defined 
around eg. (13. 141) . Then, for large r one has 



n oc N, 



Hf OC Nf 



2 3 log r - 1 



n: 



c+ - 2N,So 
6 



4 36 log r — 15 log r + 1 

6NJogr-5N,-2NfSoo 



12r2 



(5.9) 



We see that far in the UV the cascade is at work, with a number of 'bulk' D3-branes 
n = kNc — as usual k ~ Nc\og{r). The Higgsing has stopped and the source D3-branes rif 
just accounts for a constant charge. This situation is also reflected in the warp factor /i, see 
the discussion around eq.f l4.4p . 



6 Towards applications. 

This section is devoted to provide the reader with a reason, based on field-theory argu- 
ments, why this type of solutions is interesting and might be relevant for phenomenology 
and model-building in high-energy Physics. In terms of the field theory dual, the back- 
grounds constructed and discussed in this paper have a set of very peculiar features, which 
can be summarized as follows. 

First of all, notice that (besides the confinement scale related to the end of space in the 
geometry) there are three dynamical scales in the theory. One is the scale p controlled 
by the parameter c+ in eq.f l2.7p . i.e. by the insertion of a dimension-eight operator in the 
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single-site theory (which is also related to a dimension-two VEV, as explained in [19] and 
references therein). In the absence of sources, below such scale the theory is best described 
as a generalization of the one-site field theory dual to [12] , while above this scale the rotation 
procedure yields a background that is dual to a two-site quiver realizing the (baryonic branch 
of the) Klebanov-Strassler duality cascade. There are then two other scales p* and ps, 
p^, < p5, such that the function 5" has support in the range between them. There is no 
obvious relation between p and these other two scales, but for reasons of simplicity we will 
assume in this discussion that p ~ p5. Different cases can be discussed along the same lines, 
in a case-by-case way that does not add to the main physical points we want to make. 
We will focus only on the rotated solutions. 

• In the far UV, for p > ps ^ p, the theory resembles the (baryonic branch of) the KS 
cascade: the theory is flowing close to a line of fixed points, each of which is a A/" = 1 
Klebanov -Witten fixed point, a two-site quiver with gauge group SU{n) x SU{n), and 
with n increasing towards the UV (i.e. with p). The flow never approaches an actual 
fixed point at finite (and large) n because of a small imbalance between the ranks of 
the two gauge groups, but rather the flow goes up a cascade of Seiberg dualities p9] 
which continues indefinitely towards p — )■ -|-oo. See [ID]. 

• There is an intermediate range p^ < p < ps ^ p over which the function S is non-trivial. 
At the scale p and in the absence of sources, the duality cascade would stop, due to the 
Higgsing induced by the dimension-two condensate appearing and precipitating the 
theory towards the last stages of the duality cascade itself [201 [19]. On the other hand, 
because in this range the function S is non-trivial, another cascade starts, which has 
a completely different interpretation: it is a cascade of Higgsings of the gauge theory. 
Source D3s are crossed when flowing to the IR, hence realizing the ideas proposed in 

m- 

• Below the scale controlled by the value of p*, the Higgsing cascade stops, and with 
it most of the dynamical features related to Nf (up to subtleties which have been 
discussed earlier), because S vanishes. The theory now looks like a generalization of 
the single-quiver theory in [12], with the same type of dynamics. 

• At very low scales (near the end of space) of the geometry, the theory shows the 
appearance of a non-trivial gaugino condensate, and confines in the usual sense of 
producing an area law for the Wilson loop. 

Notice that one might as well take p to be very small, near the end of space: in this case the 
duality cascade would continue all the way to very small energies, in particular extending 
over the second and third of the four ranges described above. In this case, there would be a 
regime in which both Higgsing cascade and duality cascade coexist. 
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The physics taking place in the field theory at the second of the four stages (for < p < 
ps ~ p) is very peculiar, and even more peculiar and new is the fact that such a behavior 
appears only over a finite range of energies. The question is: what kind of theoretical 
models exhibit features that are qualitatively similar to the one discussed here, and what 
kind of possible physical systems can they describe? In other words: are there conceivable 
applications of such results? 

Interestingly, there is a positive answer to this question, which relates to a long-standing, 
very difficult and open phenomenological problem. A scenario that has some common ele- 
ments with our present one is that of Extended Technicolor Models (ETC) with Tumbling 
dynamics [H] . This is a very plausible dynamical explanation for the origin of fiavor physics 
(SM-fermion masses, mixing angles, CP violation, FCNC interactions, ...), in the context 
of strongly coupled extensions of the Standard Model [l2l US] (see also [S] for a sample of 
reviews on the subject). We will provide here a simple summary of what this means, besides 
referring the reader to the literature on the subject. We also want to clarify what are the 
actual similarities and the substantial differences, and hence the possible intrinsic limitations 
into attempting to use the approach of this paper in order to study tumbling ETC. 

ETC addresses the following, well-known, fundamental problem, arising in the context 
of (strongly coupled) Dynamical Electro Weak-Symmetry Breaking (DEWSB). DEWSB or 
technicolor (TC) proposes to replace the Higgs sector of the Standard Model with a new gauge 
theory with group Gtc and new (techni-fermion) fields transforming non-trivially under the 
action of both Gtc and the Standard Model gauge group. The strong dynamics associated 
with Gtc yields the formation (at the confinement scale Atc) of non-trivial condensates 
made of techni-fermions, which results in the spontaneous breaking of the SM gauge group 
at the Electro Weak scale. The SM gauge interactions themselves communicate the breaking 
to the W and Z gauge bosons, which become massive. However, in the process of removing 
the Higgs field, one also loses the Yukawa couplings, and hence one needs a mechanism that 
couples the SM fermions to the techni-quarks, in order for the quarks and leptons to acquire 
a mass below the Electro Weak scale. 

ETC provides such a dynamical mechanism. The generic ETC model works in the following 
way (an example of such a construction is discussed in detail in [15]). Start from some gauge 
theory with gauge group Getc x Gsm (where Gsm = SU{3)c x SU{2)l x U{1)y is the 
familiar SM gauge group), and a given fermionic-matter content. Assume that the dynamics 
of Getc is such that the theory is asymptotically free, but undergoes a sequence of breaking 
stages 

Getc Gi ^ G2 ^ Gtc , 

at scales Ai ^ A2 ^ A3, respectively. At scales below A3 the resulting effective field theory 
consists of the following. 

• A gauge theory with gauge group Gtc x Gsm- 
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• Two kinds of massless fermions: techni-quarks ipxc transforming non-trivially under 
Gtc X GsM, and singlets of Gtc that we denote by ipsM, that transform non-trivially 
only under Gsm- The latter are identified with the quarks and leptons of the Standard 
Model. 

• Higher-order operators originating from integrating out the heavy gauge bosons of the 
coset G ETC I Gtc- In particular, some of these are four-fermion operators coupling two 
SM fermions with two TC fermions, with the generic form 

■^tpTci^Tci'SM^^SM ■ (6.1) 
i 

We present in Figure E] a cartoon of the dynamics (in terms of the effective gauge coupling) 
of such a scenario. 




Figure 6: Cartoon depiction of the running of the gauge coupling g"^ in an example of multi- 
scale (tumbling) ETC model as a function of the renormalization scale fi. The sequential 
breaking at scales Ai > A2 > A3 > Atc means that the gauge group and field content 
changes at each breaking scale, and so does the dynamics. A model that is believed to 
exhibit such behavior is described in detail in 

Now, the resulting gauge theory is the TC theory. Ultimately, Gtc will confine, and pro- 
duce condensates that break Gsm- However, the presence of the four-fermion interactions 
means that after Electro Weak- Symmetry Breaking the quarks and leptons will become mas- 
sive. Effectively, these four-fermion operators play the same role as the Yukawa couplings 
in the Standard Model. Notice that, because they originate at different scales, there will be 
in general three families of SM fermions, and some of these operators will be suppressed as 
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others as I/A3, and so on. In general there will be a very complicated, hierarchical 
structure in the four-fermion couplings, which will translate after dimensional transmutation 
into the hierarchical structure of the masses of the SM fermions, and the mixing angles in 
the CKM mixing matrix. 

The presence of a long (at least three-stage) sequence of breaking, and of hierarchies in 
the dynamical scales is absolutely necessary on phenomenological grounds, because it offers 
the only plausible and self-contained dynamical explanation for the pattern of phenomeno- 
logical masses and mixing angles experimentally measured. On the other hand, it makes it 
extremely difficult to study, because of the strongly coupled nature of the phenomena taking 
place in these field theories. Many open questions about such dynamics, and its low-energy 
implications, require some new tools for a quantitative (and often even qualitative) analysis 
to be performed. The gravity duals we are finding resemble this scenario: a sequence of Hig- 
gsing stages taking place at several scales, over a finite energy range. It is hence conceivable 
that some of the long-standing open problems might be addressed in this context. For exam- 
ple, one would also like to understand what kind of low-energy spectrum one may observe, 
with particular reference to the presence of possible pseudo-Goldstone bosons associated 
with the breaking of accidental global symmetries. 

There is another long-standing problem (see for instance [36] and references therein). 
Besides the gauge symmetries, all of the TC and ETC models also possess very large (ap- 
proximate) global symmetries, which are spontaneously broken by the many condensates 
that form. This might yield the presence of such pseudo-Goldstones, with masses that might 
put them well below the current exclusion region, and hence render the models phenomeno- 
logically not viable. Computing masses and couplings of such particles requires dedicated 
strong-coupling calculations, and new tools are needed to perform them. 

We can hence conclude that this is a first concrete example of a model which shares some 
of the fundamental features of tumbling, and that these types of models might be used to 
characterize in field-theory terms what are the features associated with the tumbling itself. It 
must also be stressed that the very nature of these models is such that a direct comparison to 
the real world must be done with caution: effectively we have an infinite, continuum number 
of Higgsing stages, rather than a few (three) distinct and hierarchical breaking stages. In 
other words, comparing to Figure El the Higgsing cascade differs from the tumbling because 
between the scales Ai and A3 (which can be associated with ps and p^,, respectively), one 
has a continuum of breaking scales, rather than just a few steps. It hence remains to be 
understood what type of operators will replace the four-fermion operators of low-energy TC, 
and what type of phenomenological implications they have. This is somewhat analogous to 
the fact that the gravity dual of the duality cascade actually consists of an infinite continuum 
of Seiberg dualities, rather than a small number of such stages. 
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7 Conclusions 



Let us conclude this paper by recalling the main results of this work and proposing some 
topics for future investigation. The reader who wants to have a simple summary of the results 
and main original features discussed in the paper should complement this final section with 
the brief summary at the beginning of Section [6l 

The basic idea of this paper is to construct new classes of gravity backgrounds that we 
interpret in terms of the very non-trivial RG flow of a dual field theory with M = 1 super- 
symmetry, which exhibits a number of interesting properties. 

We constructed these new solutions, which generalize the KS j8] and baryonic branch [22] 
solutions, by the addition of sources into the type IIB configuration. We proposed suitable 
profiles for the sources that avoid any small radius singularity and also provide the back- 
ground with a 'logarithmically' AdS space in the large-radius regime. In this fashion, the 
backgrounds are trustable and the dual field theory associated is well defined all along the 
flow. It is therefore sensible to think that these backgrounds are characterizing a well-defined 
field theory, and hence we suggested an interpretation for the microscopic properties of such 
a field theory. We studied different aspects of the four-dimensional dual quantum field theory 
that point to the consistency of the interpretation given. 

A comment on the string side of the construction must be made. While the solutions 
with only D3 sources on the deformed conifold are fully rigorous, our profiles for D5 and D3 
sources have not been obtained from a microscopic kappa-symmetric solution to the brane 
equations of motion. In spite of this, various tests suggest that the backgrounds presented 
here are physically acceptable, both from the supergravity and from the field-theory sides. 

The main element of novelty that emerged from our study is that the dual field theory 
undergoes two cascades: a cascade of Seiberg dualities, analogous to the one in [8], and a 
cascade of Higgsings, analogous to what was suggested in p3]. In both cases, the rank of the 
gauge groups in the dual quiver field theory shrinks when going towards the IR. But this is 
the result of very different phenomena. The fact that the Higgsing cascade takes place over 
a compact, finite range of energy scales is a significant element of novelty, also considering 
that traditional field theory studies of this type of behavior are peculiarly difficult. 

There are several open questions we leave for future studies. On the more formal side, a 
well defined project would be to determine from microscopic calculations the profiles (the 
function S in the body of the paper) that reproduce the qualitative behavior we proposed 
here. This would complete the picture in a satisfactory way. Also, a side project would 
be related to a precise discussion of the existence (or not), couplings and properties of the 
Goldstone mode associated with the baryonic symmetry, a problem that we only briefly 
touched upon in this work. 

This work opens the way to several possible applications. In particular, it would be 
interesting to make more concrete the possible applications in the context of Dynamical 
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Electroweak-Symmetry Breaking, for example by calculating observables related to the dis- 
cussion in Section [HI or by adapting to these backgrounds (or suitable modifications of them) 
the tests performed in In this way, one might make some interesting progress towards 
the construction of a complete model of dynamical Electro Weak Symmetry Breaking in the 
rigorous context of gauge/string dualities, and hence move towards the direction of obtaining 
a fully computable and testable model of this type. 
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A Technical aspects of the SUSY backgrounds without 
sources 

We write in this appendix various technical aspects of the supersymmetric backgrounds 
without sources. As explained in Section [2l one can partially integrate the BPS system, and 
rewrite the functions $, h, g, k, a in terms of two new functions P and Q. This was originally 
derived in [26], and the solution one finds is: 




2 ' P cosh(2p) - Q sinh(2p) ' (P^ _ g2)p/ • 



Q 



(A.l) 



The function b{p) can be integrated exactly to give 




(A.2) 



P is given as the solution of the master equation (12. 6p . that we recall here: 




(A.3) 
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The solutions of the master equation are known only numerically, but their asymptotics can 
be written down analytically. In the UV (that is, large values of p), we had the expansion 



P=e 



4p/3 
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p + 



-20p/3\ 
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which can be plugged back into eg. flA.ip to obtain the background functions at large p, 
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The geometry in eq.f l2.4p asymptotes to the conifold after using the expansions above. In 
the IR (that is, close to the origin of the space that we take to be p = 0), we have to use 
eq.dlH]): 
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This space is free of singularities as can be checked by computing curvature invariants. 



(A.7) 



B Detailed derivations in the case of sources with pro- 
file 

Let us work out the derivation of the supergravity solutions in Section [31 for the case where 
we add sources with a profile S{r). Suppose that we start with the master equation (13 ■4p . 

-P' + Q' + 2NfS P' -Q' + 2NfS 



{P" + NfS') + {P' + NfS) 
that we will rewrite as 
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Then, we can take this last equation and integrate it twice, to get 
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We will now propose a solution in an inverse series expansion in the constant A 
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B.l The limit A ^ oo 

We describe in detail the limit of large A. This will lead to an analytic solution. As will see, 
we need to take at the same time Nf -t- to make sense of the background. First of all, we 
write the explicit expansions for the functions in the limit 

22/3 

A-^oo, Nf^O, NfX-^^u. (B.6) 
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+ 8 + ^ + ^' 

1 ^ Qtanh(2p) ^ gtanh^ ^.^^^^p) - Po cosh(2p)) + 0(A- 



cosh2p Pi cosh (2p) A Pi A 
4^-4<^o _ 2sinh^(2p) 2sinh^ 

~ A3P2P{ A4p3p;2 lAUVj^ + f^oj + 2fi3f^J + 

+ ^1^^ (^i(^o + iV/5)(2P^A + 2PoPi' + iV.PiS) - P^^Pi^+ 
+Pi'(3Po'Pi' + Pi'g' - 2PiPi'P_i - PfPii)) + 0(A-6) , 

^ kV^o / 2sinh^(2;^ I hj K , .,.,_9/2^ .R7^ 

where we have defined 



^ _/.V'^o 2sinh^(2p) 

" 2AP[ V PfP; + ^^^^^^ + 



.^e^-o 2sinh^(2p) r ^ , , ^ ^ ^ , (B.8 



8Pi'2p2 Y P2P( 



4Pi''(2Po' + g') + 4PiPi'(2P_iPi' - PoP^) + 



+ p2(4P;P:i - 3Po^) - 2A^/Pi5[3PiPo + 2PoPi'] - SNjP^S^ 



Notice that the expressions for Q{p),b{p) do not change from those in eqs. (13.31) and (13. 7p 
respectively. Now, we will consider the limit in eq.( IB.6p . We choose the value of the constant 

Kj cLS 

3k'' = 2e^A3e-^*» . (B.9) 
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Again the logic for this choice is to get a warp factor that vanishes at large p. In order to 
have a well-defined metric (where powers of A are absent) when taking the limit A — t- oo, we 
also need to take Nj — as in eq.f lB.6|) . Doing this has the following effect on the functions 
defining the full background: 



^2k 



A 



tanh(2p) 



.2g 



4 



^coth(2p)^ + 0(A°) 



4 



C(A° 



4 8 ^ ' cosh2p 
Q = Ne{2p coth(2p) - 1) + C(A"^) , b 



2p 



sinh(2p) 



3e4<fo 



(B.IO) 



and on the function h this has two interesting effects. On the one hand taking Nf — )■ 
considerably simplifies the expression for hc,hf. On the other hand the scaling NfX = 
22/3g-4/3^ fixed makes the term X~^^^hf scale like X~'^^^hc (that is as A~^). Notice also that 
in the limit of eq.( lB.6p . X~^^'^hf becomes 



x-'/% = X 



-2 



r8/3 



dx- 



S{x) 



sinh(4a;) — 4a;) 



(B.ll) 



Regarding the integrals defining he, we see that in the limit of eq. (lB.6p . we have 



he 



25/3 



X7/2 ^2^8/3 



°° , (sinh(4x) 
dx 



4x)^/^(2x coth(2x) 



sinh (2x) 



+ o{x- 



h 



KS 



A2 



+ 0{X~ 



(B.12) 

2 

Notice that we dropped terms that are suppressed like ^ in the previous expression. Finally, 
rescaling the Minkowski coordinates Xi — )■ XjA~^/^ we have a metric that is independent of 
the parameter A. Using that 

P[ ^ 8sinh^(2p) 

Pi 3(sinh(4p) - 4p) ' ^ ' ' 

the internal space metric is the deformed conifold. In this way, we have generated a new 
analytic family of solutions that depend on the function S{p). For correct choices of S{p), 
like for example S{p) ~ (tanh(2p))^" as in [3l] or those in eq. fl3.16l) . this family of solutions 



is non-singular 
eq. 



Interestingly, it is also true when we do not take the limit described in 



6]). In that case, our series expansion for P(p) and h are not truncated. If we take 
the strict Nf = case, this solution is the baryonic branch of |22]. The parameter A is the 
one moving between different VEV's for the baryon and anti-baryon operators (a parameter 
called U in [9]). For nonzero values of Nf, an interpretation of the parameter A was given in 

Eol. 
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B.2 The exact and analytic solution 



We describe here the solution obtained in the hmit A — t- oo. The warp factor is 

1/4 r°° S(x) - \ 

After rescahng Xi — )■ Xi\~^ and choosing $(oo) = 0, the full configuration will read 

ds^ = h~^'^dxl^ + h"^dsl, 

dsl = Pi (^{dp' + ]{COs + cos edpf) + ^^^M2/^) ^ gi^2 ^^^2)^ 

\2ii 4 4 

+ coth(2p) [(wi + adef + (cDi — a sin 6dp>) ] 
N^2pcoth{2p) -1 



Bo 



sinh(2p) 



cosh(2p) (sin 9d9 A dp — sin 6d6 A — sin 6dp A ui — d9 A 002 



H3 = dB2 , F3 = *6^3 , 

F5 = — (1 + *)dpH~^dt A dxi A dx2 A dxs A dp , 
H = hKs + -^i" dx ^ 



' es/s" Jp (sinh(4x) - 4x)i/3 " 

(B.15) 

Here, taking the limit for B2 is non-trivial. Indeed, when expanding it in terms of inverse 
powers of A, one gets 

B2 = X^-^^^^^^^d(p^{us + cos 9 dp)) -Bks + 0{\-') . (B.16) 

Using eg. (12. 171) . one notices that the overall factor in front of the total derivative in the 
previous expression is a constant. So the leading order term in B2 can be gauged away 
because it is exact. So one is left with B2 = —B^s + ^ {X~^)y which has a well-defined limit 
when A — )■ 00. Note that the sign difference between our result and the one of Klebanov- 
Strassler comes from a different choice of orientation. 



C Equations of motion 

In this appendix we state the equations of motion of the problem studied in this paper, that 
is type IIB supergravity with brane sources. The interested reader will find more details 
about the action from which those equations are derived in [20]. Let us first start by writing 
the Bianchi identities, which are modified due to the presence of sources. From equations 
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fl2.15p and (13. 2p . one can see that the three-form flux F3 is not closed anymore: 



Nf 

dF^ =— ^ sine de Add) A 



~,3 



Soj^ Au^ - S'dpAu 



+ 



Nf S' 



dp A 



dO Au? Au'^ + d(h A ( sin 6 u}"- A uj-^ + cos 9 dO A uj 



~,1 A ,~,3 



,~,2 



8 cosh(2p) 
We then define a four-form S4 as 

That means that we get the following Bianchi identities for the other fluxes: 

dH3 = 0, dF5 = H3AF3 + B2AE4. 

The equations of motion for the fluxes read 



1 - h 



d (e-* * H3) = F3 A F5 + e*^— ^ V0I4 A S4 , 



d (e* * Fs) = -Hs A F5 



h 



where V0I4 = dx^ A dx^ A dx^ A dx^ . We then define the following notation 
We also have that 

j w(p) A A(io-p) = - j v^Aj(*w) . 
Using these, we can write the dilaton equation of motion as 

where J is the almost-Kahler form of the internal space and is equal to 



- cosafe^^ + e^^) - sma(e'^ + e^^] 



using the conventions of equations f l2.14p - p.16p . Finally, the Einstein equation is 

I ^ ^ — ^ (ptTJ W n J-T^\ -L T? T? T1T2T3T4, I T-isources 
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p 



Figure 7: Plots of T^'?^' (^^ N^ = 2 = Nf and /ii = 5) for different profiles. The dot-dashed 
orange curve corresponds to a profile 5* decreasing like e"^^"'^-'^/^. It grows exponentially in 
the UV. The red curve is for the profile in eq. fl3.16p . and it behaves nicely. And the dashed 
blue curve corresponds to a profile S decreasing like e"*^^"'"^-*''/^. It becomes negative for large 
values of p. We take e = 1 to appreciate well the qualitative difference, but the results hold 
for any e > 0. 



where T^^^"^^^ is the energy-momentum tensor coming from the source action. It is given by: 



1 e 



2$ 



T^r^^^^ = - --^ (H,n...3 * (VoU A J)/--^ - 6^7,.S4. * (VoU A J)) - 

1 v^T^ ^^^^ ^ S4)^.,....3 (*Vol4)/^-^^ - 120(7^.(52 A S4) J (*Vol4)) . 



240 h 

(C.IO) 

From (1C.9P and (IC.lOp . one can get an equation for the Ricci tensor as: 
I ^_ F F nT2T3T4 

96 '^'^l'^2T3T4-' u 

1 (C.ll) 

- (2H^.,,,.3 * (V0I4 A J)J^^'^' - 3^?^.H4J * (V0I4 A J)) - 



1 ^#vi- ^ 

240 



e*^^ ((^2 A S4)^,,..,, (*Vol4)/^-"^ - 60(?^.(52 A ^4)^ (*Vol4)) . 



To clarify the contribution of the sources to the problem, let us express in more details 
the various components of the energy-momentum tensor due to the sources. In fiat indices. 
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the energy-momentum tensor of the sources gives: 



™rces _ Nf 2g~2h~2k-'S>/2 (2k q , 4e^^ + 6^^ (a - COsh(2p))^ \ 
T-isources ~2g—2h~^/2 ry-isources 

^sources g— 2g— 2fc— $/2^/ ^^sources /Q -j^2'l 

~ /ii/2sinh(4p) - 0-^ ' 

rpsources _ -2g-2h-2k-<S>/2 2e^^ + (q — COsh(2p)) , _ ^sources 

~ 2hy' sinh(4p) ^ -'^^ ' 

T^sources _ „-g-/i-2fc-$/2 Q — COsh(2p) , _ rynsources 

"2/^1/2" sinh(4p) ^ -^^2 • 

where rjij is the flat Minkowski metric. Notice that this energy-momentum tensor is the same 
one (up to global factors of h) as the one in the unrotated background [32]: the rotation 
does not change it. 

One expects that the tt-component of the energy-momentum tensor, T^o"o'^''^, represents the 
mass density of sources in our 'static' background. As such, it should be a positive-definite 
quantity (even if we perform a Lorentz transformation, its sign should be preserved). If we 
look at its expression from eq.f lC.12|) . we have that 



™rces Nf 2g-2h-2k-^/2 [ ^2k ^ , 4e^^ + 6^^ (o - COsh(2p)) 

2/^3/2 I smh(4p) 



(C.13) 



It is trivial that this positivity condition will hold for profiles like the one in Figure [H as 
both S and S' are positive. However, for profiles like eq.f l3.16p that we are studying in this 
paper, this condition is not trivially satisfied anymore. 

We have seen (see Figure [7]) that 5* should not decrease faster than e~^P^'^ in order for 
T^^^o'^^ to be positive everywhere. Moreover, we also know that 5* should decrease at least 
as fast as e"'^''/^ in order to preserve the KS UV asymptotics. This argument seems to imply 
that the only physical choice is actually the one we have made in eq. fl3.16p . Notice that this 
choice is the only one that naturally gives a constant density of D3 sources in the rotated 
background, in virtue of eq. (l3.1ip . 
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